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Abstract 

We discuss the current picture of the standard model’s scalar sector at 
strong coupling. We compare the pattern observed in the scalar sector in per¬ 
turbation theory up to two-loop with the nonperturbative solution obtained by 
a next-to-leading order 1/A expansion. In particular, we analyze two resonant 
Higgs scattering processes, // ^ A —> /'/' and ff ^ H ^ ZZ, WW. We 
describe the ingredients of the nonperturbative calculation, such as the tachy- 
onic regularization, the higher order 1/N intermediate renormalization, and the 
numerical methods for evaluating the graphs. 

We discuss briefly the perspectives and usefulness of extending these non¬ 
perturbative methods to other theories. 


1 Introduction 

The possibility that the electroweak symmetry is broken by a strongly interacting 
scalar sector received considerable attention in the literature. Interesting scenarios 
were proposed, such as the possibility of a Higgs boson coupled strongly to the vector 
bosons and to itself |l[], and the formation of a spectrum of bound states at a higher 
scale which would restore unitarity in scattering processes . Also phenomenological 
models were proposed for studying quantitatively the implications of strong interac¬ 
tions in the electroweak symmetry breaking sector, such as the BESS model 
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However, beyond the phenomenological models of strong interactions, an approach 
based on hrst principles was missing becanse of leak of a nonpertnrbative solution, 
and of technical difficulties in extending perturbation theory in higher-loop orders. 
Realistic calculations on a lattice of physical processes involving the Higgs sector 
are still confronted with technical limitations set, among other issues, by the size of 
the lattice. 1/N expansions in the Higgs sector were only performed at leading order, 
which is a rather poor approximation. Perturbation theory in the Higgs sector beyond 
one-loop becomes very difficult because it involves Feynman diagrams with massive 
internal lines and finite external momenta, for which already at two-loop there are no 
general analytical solutions available. 

Recently, considerable progress has been made in understanding from first princi¬ 
ples the nature of the standard Higgs sector when its coupling becomes strong. This 
is due mainly to technical advances in massive higher-loop techniques and in higher- 
order nonpertnrbative 1/N expansions. In this paper we would like to discuss the 
perturbative and nonpertnrbative aspects of this behaviour at strong coupling. 

The major question which will be addressed by future experiments at the LHC 
is how is the electroweak symmetry broken in nature. While it may or may not 
turn out to be actually broken by strong, nonpertnrbative interactions, the Higgs 
sector remains a fairly simple but not trivial model, where new perturbative and 
nonpertnrbative solutions can be tested, in view of applying them to other, possibly 
more complicated theories. 


2 Higher-order perturbation theory 


The existing calculations of leading mu radiative corrections in the standard Higgs 
sector at two-loop level are based on using the equivalence theorem in Landau gauge. 
This way radiative corrections involve only diagrams with scalars on the internal lines, 
so that the problem at hand becomes much simpler. This procedure was proposed 
for the first time in ref. , where the one-loop correction to Higgs decay into vector 
bosons was calculated in this way. 

So far at two-loop level the scalar self-energies and the main decay modes of 
a heavy Higgs boson (ft, WW, ZZ) are known, and also the high energy limit of 
vector boson scattering |^. The vector boson scattering is known completely only 
at one-loop level |^. This is due to the complexity of the diagrams involved in a 
complete two-loop treatment. By using the existing two- and three-point functions, 
some other scattering processes of phenomenological interest can be derived 0, 

For a discussion of the existing results concerning effects of enhanced electroweak 
strength in the standard model at two-loop order, see for instance ref. 0- 

The main decay modes of heavy Higgs bosons are into pairs of vector bosons and 
into top quark pairs. At leading order, these decay width are given by the following 
expressions: 
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a) b) 


Figure 1: The magnitude of the leading tuh radiative corrections to the H tf (a) 
and the H ZZ^WW (b) decays. The plots show the ratios of the decay widths 
at one-loop (solid line) and two-loop (dashed line) to the tree level decay widths as a 
function of the on-shell Higgs mass. 
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The radiative corrections of enhanced electroweak strength up to two-loop order 
are given by the following multiplicative factors [^: 
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Here A = (is the quartic coupling of the scalar sector. 
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In the above expressions, the strength of the coupling of the scalar sector is pa¬ 
rameterized by the on-shell Higgs mass mjj, which is defined by the on-shell renor¬ 
malization condition Re = m\) = 0. We plot the correction factors given 

by eqns. 2 in figure 1. One can see that for both correction factors the two-loop cor¬ 
rection becomes as large as the one-loop correction for mn about 1 TeV (1.1 TeV for 
H —> //, and 930 GeV for H WW, ZZ). Even if the higher-loop corrections are 
not yet known, this pattern suggests that the radiative corrections blow up strongly 
around mn ~ 1 TeV. Note that the scheme ambiguities associated with the two-loop 
result become substantial already for considerable smaller values of the Higgs mass 
in| . Of course, similar conclusions can be obtained from other scattering processes 
apart from the decay modes. This behaviour is in agreement with well-established 
results regarding perturbative unitarity violation in vector boson scattering at tree 
level jT^ . 

The point where the perturbative expansion blows up depends on the expansion 
parameter, and therefore on the renormalization scheme. So far the two-loop results 
mentioned above were translated into the MS []n| and the pole 0.0 renormalization 
schemes for processes at the Higgs mass energy scale, in the hope that perturbation 
theory may show better convergence properties in certain schemes. An overall con¬ 
clusion of these studies is that the MS scheme diverges somewhat sooner and has 
larger scheme uncertainty than the on-shell scheme, while the pole scheme appears 
to have slightly better convergence properties than the on-shell scheme. For a more 
rehned discussion of the scheme and gauge dependence see ref. p 

While the existing higher-order perturbative results in the scalar sector are con¬ 
sistent with a strong blow-up of radiative corrections at about rriH ~ 1 TeV, this 
fact appears somehow puzzling if one considers that the quartic coupling of the scalar 
sector, A = is only of order .4 for uih ~ 1 TeV. 

From the perspective of the nonperturbative solution to be discussed in the fol¬ 
lowing section, the reason for this behaviour is the Higgs mass saturation effect. Due 
to the dynamics of the scalar sector, the on-shell Higgs mass is not a good param¬ 
eterization of the quartic coupling for values larger than mn ~ 1 TeV. However, as 
will be shown in the following section, at the fundamental level, the scalar sector is 
perfectly well dehned at higher values of the quartic coupling. Well-behaved, unitary 
solutions can be obtained for various processes by using nonperturbative methods. 

At a more speculative level, the fact that the pole scheme seems to show somewhat 
better convergence properties than the on-shell scheme may have its origin in another 
nonperturbative property of the Higgs self-energy. This type of nonperturbative effect 
was discussed in ref. [113 in the context of a theory of vector bosons coupled to a 


large number of light fermions. The idea is that due to the structure of the self-energy 
function the on-shell renormalization condition may have no solution for a range of 
the theory’s coupling, while the pole renormalization condition may have a solution. 
This results in the on-shell mass being suitable for parameterizing the theory only 
for a more limited range of the physical coupling than the pole mass. Should such a 
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nonperturbative effect be present in the Higgs sector, this conld explain the pattern 
observed in pertnrbation theory. In ref. this was investigated by using the 


nonperturbative 1/N expansion, and it was shown that this effect is not present in 
the Higgs sector of the standard model at LO. It is still an open question if it appears 
at NLO. 


3 Nonperturbative 1/N solution 

The 1/N expansion aims at a nonperturbative solution in order to avoid the problems 
of perturbation theory at large coupling. Perturbation theory ceases to be a satisfac¬ 
tory solution when radiative correction blow up already at lower loop orders and the 
renormalization scheme ambiguity is so large that the result becomes unreliable. 

The 1/N approach is free of these problems because all radiative corrections of 
all loop orders are explicitly summed up. The idea is to treat the Higgs sector as an 
0(V)-symmetric sigma model, where N = A for the standard model, and to expand 
in 1/V instead of the quartic coupling. The solution is then valid independently of 
the strength of the coupling — the quality of the approximation depends on the value 
of N. Also it is completely free of renormalization scheme ambiguities. One can work 
in any intermediate renormalization scheme and still obtain the same result. 

Another interesting nonperturbative feature of the 1/N expansion of the sigma 
model is the hniteness of wave function renormalization constants. This is a property 
of the exact nonperturbative solution, and was checked at next-to-leading order in 
1/N. However, the renormalization of the coupling constants is ultraviolet divergent. 

As an extra bonus, the 1/N solution provides naturally a consistent treatment of 
resonant scattering amplitudes. This was known as a long-standing issue in perturba¬ 
tion theory. The essence of the problem is that around a resonance one has to perform 
a Dyson summation which in perturbation theory at any hnite order introduces in¬ 
complete higher order contributions which are unphysical. In gauge theories this leads 
to gauge dependent results. Meanwhile a solution to this problem was found PHH, 
which is based on a Laurent expansion around the physical pole. Within perturbation 
theory, this solves the problem in a fundamental way, and applies consistently to all 
orders in perturbation theory. Another approach which was proposed is to use gauge 
invariant pinch technique self-energies in the Dyson summation. For phenomenolog¬ 
ical purposes only, other approaches to treat resonant amplitudes were proposed in 
the literature, which amount to special resummations of low oder results. Some of 
them may be easier to use in certain phenomenological applications. At the same 
time they are less fundamental theoretically — some of them still contain unphysical 
higher order terms, or apply only to tree level or one-loop calculations. 

The 1/N solution is automatically free of these problems of perturbation theory 
because it is an all-order solution in the loop expansion. 

However, the 1/N solution still has some residual ambiguity which can be related 
to the triviality problem and to possible nondecoupling effects from a hidden heavy 
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sector. Technically, this appears in the tachyonic regnlarization. Pertnrbation theory 
is used at an intermediary stage in the usual 1/N treatment, but perturbation theory 
does not determine the solution uniquely. This is the physical origin of the ambiguity 
entailed in the tachyonic regularization. 

3.1 1/N combinatorial rearrangement and diagrammatics 

We start with the usual Lagrangean of a 0(W)-symmetric sigma model: 

A = i4'S>o9''4o - f ■S-S - ■ ■ ■ A«) (3) 

From this Lagrangean one can in principle derive directly a perturbative expansion 
for Green functions, and classify the Feynman diagrams according to their order in 
1/N. However, beyond leading order the combinatorics becomes very complicated. 
In order to perform explicit calculations beyond leading order it is useful to perform 
a rearrangement of perturbation theory. A useful trick for doing this was proposed 
in ref. It consists of adding a nondynamical term to the Lagrangean: 

= -d^^od - -Xo^l + ^Xo - 

This involves the introduction of an unphysical auxiliary field y. As one can 
see, the equation of motion for x is just an equation of constraint, and therefore 
the physical spectrum and the dynamics of the model remain unchanged. The effect 
of this trick is that the Feynman rules are changed. Namely, the quartic couplings 
are eliminated. The only vertices left are trilinear, and involve one y field and two 
physical scalars. This simplifies enormously the combinatorics of Feynman graphs in 
higher orders. 

We note that some other rearrangement schemes with different properties were 
discussed for the 0(W)-symmetric sigma model in ref. [^. To our best knowledge, 
these schemes were not applied so far in actual calculations. 

In the following we describe the counterterm structure which is used for performing 
renormalization at NLO in the 1/N expansion. A somewhat different approach is 
presented in ref. [^, but the hnal results are the same. 

Renormalization is performed in principle order by order in perturbation theory. 
However, for performing actual calculations of higher order in 1/N, it is of advantage 
to group all counterterms of various loop orders and which are of the same 1/N order 
into the same 1/N counterterm. We define the 1/N counterterms as follows: 
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Here we already used the fact that the tadpole and wave function renormalization 
counterterms do not receive contributions at leading order in IjN. We also note that 
although two tadpole counterterms are present, bt and they are related through 
the gap equation |^, for instance by requesting that the leading order ground state 
condition y = 0 be preserved in higher orders, where y is the vacuum expectation 
value of the y field in the spontaneously broken phase. 

At this point it is useful to note that since the two Lagrangeans of eqns. 3 and 4 are 
equivalent, a linear combination will also describe the same physics. This observation 
can be exploited for performing BPHZ renormalization in a more elegant way, as will 
be explained in the following. Beyond leading order in 1/A^ it is advantageous to work 
with a linear combination of the potential parts of Lagrangeans L\ and £ 2 - Keeping 
only the contributions relevant for next-to-leading order calculations, we consider in 
fact the following Lagrangean: 
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Here K is in principle a completely arbitrary constant. We have the freedom to 
choose it so that actual calculations are more convenient. We will consider K to 
be of order 1 in the 1/N expansion. Thus the potential part of £2 is regarded as an 
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Figure 2: 


The leading order bubble diagram. 


0{1/N) counterterm. We will choose the actual value of K so that the renormalization 
procedure is more transparent at NLO in the 1/N expansion. 

The Feynman rules can be read out directly from the above expression. Within 
this diagrammatical rearrangement of the sigma model, counting powers of 1/N in 
multiloop diagrams is straightforward: closed Goldstone loops contribute a factor iV, 
while XX propagators give a 1/N factor, and mixed propagators contribute 1/\/N. 
At the same time, the absence of quartic couplings at tree level reduces considerably 
the number of possible topologies. 

For these reasons, for a given process and a given order in 1/N, it is easy to 
write down the Feynman graphs of all loop orders. As it will become clear from 
the discussion of two- and three-point functions, there is always a hnite number of 
multiloop topologies, where one can only insert chains of one-loop bubbles in the xx 
and xo" propagators without increasing the 1/N order of the graph. 

We emphasize that this combinatorial rearrangement of the sigma model is quite 
crucial. It is possible to calculate explicitly nonperturbative processes in the Higgs 
sector precisely because the combinatorial rearrangement enables one to write down 
explicitly and in a manageable way the diagrams of all loop orders, without truncating 
the perturbative expansion. 


3.2 Tachyonic regularization 


It is straightforward to derive the two-point functions of the theory at leading order 
in 1/N. This was done for instance in ref. [|^. The only diagram involved is the 
one-loop bubble diagram shown in hgure 2. One hnds the following leading order 
propagators: 
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Here is the ultraviolet finite part of the self-energy diagram 

of fig. 1, and /i is the subtraction scale. 

In the expressions above, apart from the expected Higgs pole, one notices the 
presence of a tachyonic pole. It appears at an energy s = —which is given by the 
following transcendental equation; 
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The tachyon scale A* differs from the Landau scale Al = by a shift of order 

vVAl 

The leading order tachyon is a well-known difficulty of the 1/N nonperturbative 
treatment of the sigma model. From the technical point of view, it induces causality 
violating effects in the theory. As long as one is concerned only with the leading 
order, one can try to make sense of the result by limiting its validity to an energy 
range considerably smaller than the tachyon scale. However, there is no such easy 
way out for calculations beyond leading order because the tachyon appears then in 
loops. 

One way to circumvent this is by making the assumption that the tachyon indicates 
the triviality of the theory. This sets a limit for the validity of the 1/N result at high 
energy. At some energy scale new physics sets in. This scale is presumably of the order 
of the tachyon scale, but not necessarily equal or lower. Then an obvious treatment is 
to introduce a cutoff in the loop integrations . This is then interpreted as a model 
of nondecoupling effects from an unknown heavy sector. However, in this approach 
the momentum cutoff has to be lower than the tachyon scale, which is necessary for 
computational purposes only and is not motivated physically. Also a loop momentum 
cutoff spoils the gauge invariance of the gauged model. It also introduces quadratic 
dependencies on the cutoff scale and these are known in effective theories not to be 
directly related to heavy mass effects — actually counterexamples were found in the 
literature in two-loop calculations . 

We use a different treatment of the tachyonic pole which is more convenient 
for higher order calculations in the 1/N expansion. We subtract the tachyon min¬ 
imally at its pole, which means using the following propagators instead of those of 
eqns. 7: 
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is the residuum of the tachyonic pole. 

The justihcation of the tachyonic regularization introduced above is the following. 
Green functions — such as the two-point functions above — are calculated in the 
1/N expansion starting with the perturbative expansion in the coupling constant 
A of the coefficients of the 1/N expansion. Then, all Feynman graphs of all loop 
orders which contribute to a given order in 1/N are calculated explicitly and summed 
up. Since the 1/N coefficient is only known as a power series in A to start with, 
it will be determined by its perturbative expansion up to a function of A which 
vanishes identically in perturbation theory, of the type Since the residuum of 

the tachyon is such a function whose perturbative expansion vanishes, its presence 
cannot be taken seriously as a prediction of the theory and as an indication that the 
theory is ill-dehned. While the A0^ theory is widely believed to be trivial, the tachyon 
in the 1/N expansion is certainly not a rigorous proof thereof. 

When we use the minimal tachyon subtraction scheme of eqns. 10, we effectively 
use the freedom to add a function which vanishes in perturbation theory for restoring 
the causality of the theory. In this sense, our treatment is independent of whether 
one regards the theory as being an effective one or not. If one takes the view that 
the theory is trivial and wants to include nondecoupling effects from a heavy sector, 
such effects can be superimposed over the whole calculation. Compared to a naive 
cutoff, this tachyonic regularization does not require the heavy sector to be strictly 
under the tachyon scale. 

3.3 Nonperturbative two- and three-point functions at NLO 

Beyond leading order in 1/iV, actual calculations have to be performed numerically 
because in general the multiloop diagrams involved are not manageable analytically. 
This brings about some technical complications related to the treatment of ultraviolet 
divergencies in conjunction with numerical integration. 

A useful observation is that the final result in the 1/N expansion is free of any 
renormalization scheme ambiguity. This is because it is exact at all orders in the 
coupling constant. This leaves us the freedom of working in any intermediate renor¬ 
malization scheme at our convenience, since the final result is independent of that. 
This can be best exploited for simplifying to some extent the numerical work. 

The graphs needed for the calculation of all two- and three-point functions of 
the theory at next-to-leading order in 1/N are shown in £g. 3. In addition to the 
two- and three-point graphs, there is also one tadpole graph which is needed for the 
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Figure 3: Infinite sums of multiloop Feynman diagrams which contribute in next-to- 
leading order in 1/N to the two- and three-point functions of the 0{N) sigma model. 
The blob on propagators denotes the summed-up leading order propagators. Note that 
the Tin propagator at leading order in 1/N is a free propagator. One of the graphs 
above is shown in expanded form in fig. 4- 
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Figure 4: Multiloop diagrams with three-loop topology which contribute to the XX 
propagator in next-to-leading order. 


determination of the tadpole counterterm St^. Each graph is in fact a sum of multiloop 
Feynman graphs which are all of the same order in 1/N, and of various orders in the 
coupling constant A. This is shown explicitly for one particular self-energy graph in 
fig- 4. 

The graphs involved, Ai, Ci, Di, Ei, Fi and T, are in general ultraviolet di¬ 
vergent. Since our strategy is to calculate them numerically, we hrst subtract the 
divergencies and sub divergencies of these graphs. An inspection of the Feynman dia¬ 
grams which compose each graph of hgure 3 reveals that the ultraviolet divergencies 
are polynomial, just as if they were usual Feynman diagrams, and in spite of the 
inhnite number of loops involved. We dehne in hgures 5, 6 and 7 a set of ultravio¬ 
let subtracted graphs, Aj, i?j, Q, £)j, Ej, and Fj. They are hnite, and thus can be 
calculated by direct numerical integration. 

The numerical evaluation of the ultraviolet hnite, subtracted graphs is done by 
using a numerical method for the calculation of massive three-loop Feynman dia¬ 
grams . This method reduces all subtracted graphs to a two-dimensional integral 
representation. After an appropriate rotation of the integration path in the complex 
plane, these two-dimensional integrals can be evaluated numerically ||23|| . 

The subtracted graphs are used further for calculating physical amplitudes. Here 
we consider the scattering processes // —> —> f'f' and // —> if —^ ZZ, WW . Phe¬ 

nomenologically they are important as a Higgs production mechanism at a possible 
muon collider. Also the heavy Higgs effects in these scattering processes are related 
to those in the gluon fusion process, which is the main Higgs production mechanism 
at the LHC. 

To start with, we consider the matrix element for the fermion scattering process 
//—>■//—>■ f f at leading order in the fermion mass. Then, the correction to the 
Hf f vertex is given simply by the ratio of the wave function renormalization factors 
Zo-Z^TT, because true vertex diagrams are of higher order in the fermion mass. The 
calculation reduces essentially to evaluating the Higgs propagator. Up to the overall 
factor from the tree level Yukawa couplings, the amplitude is given by the following 
expression, were we included the relevant counterterms in addition to the pure 1/N 
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Figure 5: Definition of the subtraeted xx self-energy graphs. The symbol d indicates 
the differentiation with respect to the external momentum of the box. 
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Figure 6: Definition of the subtracted aa, tttt self-energy graphs. 

multiloop graphs of fig. 3: 
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( 12 ) 


One can easily see in the above expression that one could have done without a 
wave function renormalization for the unphysical y held, and also without the K 
counterterm which was introduced in eq. 6 only for convenience purposes. It can 
be seen easily that these terms cancel ont trivially in the expression above. In fact, 
since there is no Higgs external leg in the process considered, it is also unnecessary 
to introdnce a wave fnnction renormalization for the a held. 

As we discussed in section 3.1, these connterterms are only introdnced for conve¬ 
nience. If one cancels ont these spurious terms in the expression above, one is left 
with a snm of multiloop diagrams which are individually ultraviolet divergent. Since 
the whole expression is a physical qnantity, the ultraviolet divergencies cancel out 
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Figure 7: Definition of the subtraeted vertex graphs. 


among the multiloop diagrams. However, the actual cancellation pattern is not very 
transparent due to the complexity of the diagrams. Because we need to calculate 
the 1/N graphs numerically, we need the expressions to be ultraviolet convergent. 
At the same time it is very complicated to extract the ultraviolet divergencies and 
sub divergencies from the graphs of £g. 3 as 1/e poles, as it is done in usual Feynman 
diagrams. In more complicated processes, such as three- and four-point processes, 
the cancellation is even more involved. 

The K and 5Z^ counterterms serve as vehicles of the ultraviolet cancellations 
among the multiloop diagrams. For this purpose we assign to these two counterterms 
the following expressions: 

K = -\bM 

~ < 71 ( 0 )-|-( 7 * 2 ( 0 )-|-2i?2(0)Q;o(h^ 0) , (13) 

and also we note the identity: 

5t^ = BfiQ)-BfiQ) . (14) 

Then the actual multiloop 1/N graphs from a, f3 and 7 combine with the coun¬ 
terterms and give precisely the subtracted multiloop graphs Aj, Bi and Ci dehned in 
hgures 5 and 6 : 
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/3(s) - {SZ, - SZ^) 

Actually the finite term in this expression simply means that one has to 

subtract the momentum derivative of diagram D{s) from diagram -Bi(s), rather than 
the derivative of -Bi(s), as defined in figure 6. The finite contribution which is 

left in eq. 15 simply reminds that for specifying the strength of the coupling of the 
theory, a mass scale needs to be given along with the value of A. A shift in can 

be absorbed into a shift in the subtraction point /r. As such, can be shifted to 

zero. 

Along the same lines, the following expression is obtained for the amplitude of the 
// H ^ ZZ, WW scattering process: 


/i(s) = 


m?(s) 


d(s) + 5 A^*"- 1 + 27 (s) + 


m?(s) 


Mww 

/i(s) 

/2(s) 


rn?{s) 


1-^/2 


\/Nvs — m?{s) l — 
m^(s) 


d(s) + 27 ( 5 ) + 




/3(s) 


s — m?{s) 


{5Z, 


m^(s) 


d(s) + 7 (s) - 0 (s) 




fi{s) 



All quantities involved in this expression, a, /3, 7 , fj, and <p, are sums of the 
subtracted graphs defined in figures 4, 5 and 6 . They can be calculated directly by 
numerical integration p3|, |T^ . 


3.4 The saturation effect 

The shape of the Higgs resonance can be obtained nonperturbatively in the quartic 
Higgs coupling at next-to-leading order in the 1/N expansion by evaluating numer¬ 
ically the expressions for M.fj and M.ww given in the previous section. These two 
scattering processes are the main production and decay modes for the Higgs boson 
at a possible muon collider. Also these processes are related to Higgs production by 
gluon fusion |^. This is the main production mechanism at hadron colliders such as 
the LHC. 

We give in figure 8 the resulting line shapes of the Higgs resonance. One feature of 
these line shapes is that they agree remarkably well with the perturbative results for 
low couplings. The next-to-leading order 1/N results are very close to the two-loop 
perturbation theory line shapes for mu up to about 800-900 GeV. The agreement 
confirms the consistency of the 1/N approach in higher orders and establishes that 
the next-to-leading oder is an excellent approximation. 
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Vs (TeV ) 

Figure 8: The line shape of the Higgs resonance in the scattering processes f f —> 
H —> /'/' and f f H ^ ZZ, WW. We marked the position of the maxima of the 
resonances (solid line for the 1/N result and dashed line for the perturbative result at 
two-loop). 
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At higher quartic coupling a saturation effect sets in. The maximum of the res¬ 
onance does not shift towards higher energy, the resonance only becomes wider. To 
make this effect clearer, we extract some mass (Mpeafc) and width {Tpeak) variables 
from the position and the height of the resonance in fermion scattering. The dehnition 
of the variables Mpeak and Tpeak is the following. We determined numerically the posi¬ 
tion and the height of the maximum of the resonance. Then, Mpeak and Tpeak are the 
mass and width of a Breit-Wigner resonance which has the same height and position 
of the peak. Of course, the actual Higgs lineshapes are not exactly of Breit-Wigner 
type. However, Mpeak and Tpeak describe reasonably well the main features of the 
resonance. We compare in hgure 9 the Mpeak - Tpeak relation with the perturbative 
result. Of course one is free to choose any other parameterization of the resonance, 
but the Mpeak and Tpeak variables which we use here are sufficient for comparing with 
perturbation theory. 

As one can see in hgure 8, the saturation effect is present in the ff —>■ H —>■ 
ZZ, WW scattering process as well, at a comparable energy. The precise maximum 
position of the peak is process dependent because the resonance shape is deformed 
by the energy dependence of different contributions, such as the vertex corrections in 
this case. A universal way of parameterizing the saturation effect would be by using 
the pole mass and width of the Higgs particle. Extracting this in the next-to-leading 
1/N approach is numerically more difficult because the subtracted 1/N graphs need 
to be continued into the second Riemann sheet for solving the pole equation. This is 
not available yet. 

The saturation effect provides more insight in the way perturbation theory breaks 
down and radiative corrections blow up in the Higgs sector. The dramatical failure 
of perturbation theory at around 1 TeV is well-established by tree level unitarity 
violations |^. Higher order radiative corrections blow up at a similar scale, as was 
discussed in section 2. Less violent problems can show up in the form of considerable 
scheme ambiguities of the results. At the same time however, it is puzzling that 
at the 1 TeV scale the quartic coupling is numerically not exceedingly large yet. 
The quartic coupling only becomes of the order of unity when the tree level on-shell 
Higgs mass is of the order of 1.5 TeV. Naively this is where one would expect heavy 
problems for perturbation theory to set in. From the perspective of the saturation 
effect, the explanation is clear now. The mass of the Higgs boson is not a good 
expansion parameter beyond the saturation point. In the saturation region the mass 
does not increase when the coupling is enhanced. In fact, the width of the resonance 
as a function of the mass is a double valued function. The well-known perturbation 
theory results break down at about 1 TeV because of the use of a bad expansion 
parameter. For that region the width would be a more appropriate parameterization 
of the coupling than the mass. 
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Figure 9: The saturation effect in the f f ^ H ^ f f scattering. The mass and 
width variables Mp^ak and Fpeafc are related to the position and the height of the Higgs 
resonance as explained in the text. 





4 Nonperturbative solutions of other theories: 1/N 
QED 

In this section we would like to comment briefly on the perspectives of extending this 
kind of nonperturbative methods to solving other held theories of physical interest. 

We already mentioned that the applicability of the nonperturbative 1/N expan¬ 
sion, especially in higher orders, depends crucially on the particular theory under 
consideration. One needs to arrange the perturbative expansion of the theory in such 
a way that the various loop contributions can be sorted out by powers of 1/iV in a 
manageable way, so that the graphs of all loop-orders can be explicitly calculated and 
summed up for a given order of 1/N. 

Such a theory where the rearrangement of perturbation theory is straightforward 
is the ordinary QED. QED can be seen as an example where the 1/N expansion works 
poorly because the value of N is too small. The existing calculations in perturbative 
QED are rather advanced, the coupling constant is small, and therefore the 1/N 
treatment of QED can hardly compete with the perturbative treatment. Here we 
will briefly discuss a simple leading order 1/N QED calculation for comparing with 
perturbation theory, and see in which cases a 1/N treatment can be superior to 
perturbation theory. 

QED can be organized as a 1/N expansion by introducing N species of electrons 
and at the same time dividing the gauge coupling by \fN. Then one sees immediately 
that the counting of powers of 1/A^ proceeds similarly as in the case of the sigma 
model. Ordinary QED is recovered in the limit N = 1. Already at this point one can 
expect the convergence of the 1/N expansion to be poor because of the value of the 
expansion parameter. 

We show in figure 10 the Feynman diagrams which contribute to the anomalous 
magnetic moment of the electron in leading order of 1/A^. They are the same diagrams 
discussed in ref. in the context of the large order behaviour of QED. We refer to 


this work for details on the calculation of the contributions to the anomalous magnetic 
moment due to these diagrams. The equivalent of the tachyonic regularization in the 
0(A^)-symmetric sigma model is the subtraction in the loop momentum integration 
of the diagrams of figure 10 of a term of the type which vanishes in perturbation 
theory, as discussed in ref. p3. Note that due to the value of the coupling constant. 


in QED the Landau pole is at a very high energy scale. 

The perturbative results for the anomalous magnetic moment of the electron up 
to four-loop are 


^1 — loop - 

.5 

^2—loop 

-.328478965 

^3—loop 

1.181241 

^4—loop — 

-1.557 , 


( 17 ) 
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Figure 10: Diagrams which contribute to the anomalous magnetic moment of the 
electron in leading order of the 1/N expansion. 

and their contribution to the anomalous magnetic moment is T^n an{o:/ tt)^. 

By numerical integration one can calculate the LO contribution in the 1/N ex¬ 
pansion to the anomalous magnetic moment of the electron. The result is = 

.001161494. The hrst few terms in the loop expansion of this result are: 

= .5 

= .015687421 
= .002558524 
= .000876865 
= .000470905 
= .000344687 
= .000318119 
= .000352804 
= .000455370 
= .000668824 
= .001099482 
= .001997590 

= .003971440 , (18) 


LO l/N 
^1—loop 
LO 1/N 
^2-loop 
LO 1/N 
^S—loop 
LO 1/N 
^4—loop 
LO 1/N 
^5—loop 
LO 1/N 
^d—loop 
LO 1/N 
^7—loop 
LO 1/N 
^S—loop 
LO 1/N 
^9—loop 
LO 1/N 
'^10-loop 
LO 1/N 
^11-loop 
LO 1/N 
'^12-loop 
LO 1/N 
'^13-loop 
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such that their contribution to the anomalous magnetic moment is 
The numbers above agree with the results given in ref. . 


,LO 




vr 


By comparing the LO 1/N result with perturbation theory, one sees that this 
is only marginally different from the one-loop perturbative result. Comparing the 
perturbative expansion of the LO 1/N result (eqns. 18) and the perturbation the¬ 
ory result (eqns. 17), it can be seen that the LO 1/N Feynman diagrams are not 
numerically the main contribution from higher-loop orders in perturbation theory. 

This shows that the usefulness of the 1/N expansion versus perturbation theory 
depends crucially on the value of N and of the coupling constant. If N is large 
enough, the other diagrams appearing in higher-order perturbation theory besides 
those included in the 1/N expansion will be suppressed by powers of 1/N. This way 
the 1/N expansion can converge faster than perturbation theory. If the coupling 
constant is large, perturbation theory starts to diverge already at low orders and 
the scheme ambiguities may prevent one to obtain the desired accuracy. Then the 
1/N expansion is expected to be a better alternative because it is free of scheme 
ambiguities and the solution is valid for strong coupling as well. 

Therefore it is not surprising that the 1/N expansion works so well in heavy Higgs 
physics, where = 4 and A is of order 1. 

Extending these computational tools in the case of nonabelian gauge theories is 
very difficult because of the presence of trilinear and quartic couplings of the gauge 
bosons. In the case of QCD the topological structure of the graphs which appear in 
the 1/Nf. expansion is known [^, but so far they could not be calculated explicitly in 
four dimensions. The well-known solution of QCD in two dimensions was possible 
because in this case the quartic and trilinear couplings are absent, which reduces the 
types of topologies of Feynman graphs appearing in higher loop orders. 


5 Conclusions 

We investigated in some detail the scalar sector of the standard model at strong 
coupling. For doing this we used both perturbation theory up to two-loop order 
and a nonperturbative treatment within the next-to-leading order 1/N expansion. 
With these approaches we treated the main heavy Higgs decay modes as well as two 
scattering processes where the Higgs line shape can be observed, // —> FT —>• f'f' 
and ff ^ H ^ ZZ, WW. These two scattering processes are the main s channel 
production modes of the Higgs boson at muon colliders, and can also be related to 
the heavy Higgs effects in the gluon fusion process, which is the dominant production 
mechanism at the LHC. 

The results show in all cases a very good agreement between perturbation theory 
and the 1/N expansion up to 800-900 GeV. They confirm the existence of a Higgs 
mass saturation effect in both scattering processes analyzed, which seems to be a gen¬ 
eral feature of resonant Higgs processes. The nonperturbative mass saturation value 
is just under 1 TeV for both processes. As a preliminary study based on perturbation 
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theory has shown, it is expected that the nonperturbative mass saturation effect will 
play an important role in the experimental strategy for heavy Higgs searches at the 
LHC. 

A comparison of the nonperturbative behaviour at strong coupling to the pattern 
observed in perturbation theory in higher loop orders suggests that radiative correc¬ 
tions in the scalar sector blow up very strongly at about 1 TeV mainly because of 
the nonperturbative mass saturation effect, rather than because of a genuinely strong 
coupling. Due to this dynamical effect the mass of the Higgs boson, although widely 
used as a parameterization of the coupling in phenomenological studies so far, is not 
an appropriate parameter in the saturation region. 

On the theoretical side, we have shown that in the case of the 0(iV)-symmetric 
sigma model, a reliable nonperturbative solution can be obtained by using a higher 
order 1/N expansion, which is free of renormalization scheme ambiguities and which 
is valid at strong coupling as well. We discussed the tachyonic regularization which 
we introduced for calculating higher orders in the 1/N expansion. We also discussed 
the renormalization within the 1/N expansion in higher orders. 

Finally, we discussed briefly the applicability of such nonperturbative methods to 
other theories of physical interest. 
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